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I. 



INTRODUCTION 



In recent years several mathematical models have been 
developed to predict the atmospheric dispersion of pollutants 
emitted from aircraft-related activities at and around airports. 
These models have used the steady state Gaussian plume formu- 
lation. The Gaussian formulation is used because it is adapt- 
able to distances and pollutant travel times associated with 
airports. An early contract sponsored by the U.S. Environmen- 
tal Protection Agency (EPA) resulted in a model being developed 
by the Northern Research and Engineering Corporation (Ref. 1). 
Thi3 model was later modified by GSCMET, Inc. (Ref. 2), and 
dealt specifically with civilian airport operations. A more 
recent model has been developed by Argonne National Laboratory 
(ANL) for the USA? and was termed the Air Quality Assessment 
Model for Air Force Operations (AQAM) (Ref. 3)* This computer 
model was based upon an earlier TRW model, the Air Quality 
Display Model (Ref. k) • 

Each of the models utilizes a method for solution of dif- 
fusion equations assuming Gaussian dispersion in both the hori- 
zontal and vertical directions. Gaussian formulation in air 
quality model calculations requires meteorological inputs 
including stability of the atmosphere, mixing layer height, 
and wind direction and speed. Detailed pollution source data 
are also required. The resultant models consisted of emission 
and dispersion programs. AQAM included three major parts, a 
Source Inventory model which yields annual emission at an 
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activity by source, a Short Term dispersion model which performs 
hourly-averaged calculations using input dispersion parameters 
and a Long Term dispersion model. The models predict average 
steady-state concentrations during the specified time interval 
over a specified grid surrounding the airport. 

Model verifications have to be conducted to test the algo- 
rithms and plume dispersion equations. Initial efforts to 
validate AQAM were begun by the Air Force at Williams AFB, 
Arizona. Williams AFB was chosen because it was a high traffic- 
volume, military airfield where accurate statistics would be 
available. These statistics included aircraft type, mix, and 
activity schedules from which emissions input data could be 
calculated (Ref. 5)* The objectives of the validation effort 
were three-fold: 

1. Collect a data base of airport-related air quality 
measurements to evaluate the Air Force AQAM model. 

2. Determine the impact (if any) of airport-related 
activity on local (5 km radius) air quality. 

3. Conduct a series of special studies to provide infor- 
mation on horizontal and vertical dispersion to 
supplement any model revision by ANL (Ref. 6). 

The Navy became interested in the Argonne model capabili- 
ties relative to Naval Air operations. Under sponsorship of 
the Naval Air Propulsion Center (NAPC) Trenton, N.J., the Naval 
Postgraduate School (NPS), Monterey, Ca., obtained copies of 
both the Source Inventory and the Short Term models of AQAM 
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for evaluation and adaptation to Navy operations. Upon com- 
pletion of modifications, a validation effort similar to the 
one at Williams APB was planned at NAS Miramar, California. 

The Source Inventory Program, as originally received from 
ANL, computes annual emissions of three types of sources: 
aircraft, airbase (non-aircraft) and environment (off-airbase). 
Bach of these types is further reduced by geometric configura- 
tion to either a point, line or area source. Data are input 
to the Source Inventory program relative to the type and size 
of source, location of the emission plume in three-dimensional 
space and the mass emission rate of each pollutant emitted by 
the source. The model input is often comprehensive and volum- 
inous, leaving a great margin for possible error. The program 
calculates annual emissions and provides a qualitative ranking 
of the contributions to the ambient air pollution of any indi- 
vidual source. It also prepares a data bank containing source 
characteristics, annual emission rates and temporal distribu- 
tion activity for utilization by the Short Term program. 

The Short Term program receives the above compiled annual 
results and calculates the dispersion of generated pollutants 
over a specified receptor grid during a given hour, day and 
month utilizing average meteorological data input for that 
hour (Ref. 7). For point and area sources this is accomplished 
by using initial source dimensions and meteorological stability 
criteria to project a pseudo-upwind point source. Line sources 
are generated along the route of travel of the source vehicles. 



The Short Terra model utilizes a line dispersion theory devel- 
oped by ANL. The line of finite cross-section is segmented 
into shorter lines, or "puffs", which are then dispersed from 
pseudo-upwind line sources in much the same manner as point 
and area sources (Ref. 3»3)» 

Principal modifications to AQAM were required by the Navy 
due to differences in flight operations between the Navy and 
Air Force. Subroutines were added to AQAM to account for 
Visual Flight Rule (VFR) approaches including aircraft entry 
break above the runway. Navy touch-and-go cycles, field carrier 
landing practices (FCLP), takeoff delays, and hot refueling 
(refueling of aircraft while engines are operating). Also 
AQAM was expanded to handle helicopter operations. It should 
be noted that modifications were only made to subroutines in- 
volving aircraft sources. Airbase and environ source data 
remain relatively consistent from base to base whether Navy, 

Air Force or civilian. The Short Term portion of AQAM was 
modified to calculate dispersion of pollutants over lj.12 grid 
receptors rather than the Air Force’s 312 receptors. This was 
done so that a larger off airbase area could be included in 
the analysis. Finally, Navy aircraft engines and fuel types 
are often different than those of the Air Force and, conse- 
quently, aircraft performance data and emissions data had to 
be input to reflect the changes. A plot routine was also 
incorporated into AQAM so that predicted pollutant distribution 
patterns could be more readily observed (Ref. 9). 
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The aforementioned model verification performed by the 
Air Force at Williams AFB involved 1 3 months of continuous 
air monitoring during the period June 1976 through June 1977. 
Air quality data were collected at five ground stations and 
meteorological data were taken routinely at the base weather 
station. Aircraft operations data and airbase and environ 
source information were then input to A QAM and predicted 
values of pollutant concentrations were compared with observed, 
or measured data from the monitoring stations. Preliminary 
results have indicated that a reasonable correlation exists 
between predicted and observed hourly pollutant levels (Ref. 
10 ). 

The Air Force effort included a wide range of meteoro- 
logical conditions collected over a long period of time. It 
was decided to concentrate the Navy validation effort on a 
specific meteorological "window" which would be reasonably 
stable for several days and which would occur when a large 
amount of aircraft activity occurred. The latter was neces- 
sary in order to minimize the problem inherent with high back- 
ground pollution levels. Specifically, it was desired to 
perform the validation effort at NAS Miramar, CA and to obtain 
more detailed data relative to (1 ) aircraft taxi and refueling 
operations, (2) hourly aircraft flight activity, and (3) 
meteorology. 

Once the Navy modifications were completed and input data 
were obtained for NAS Miramar, it was necessary to determine 
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the sensitivity of the model predictions to the input meteo- 
rological and operational conditions and to certain dispersion 
model parameters. An initial model sensitivity study was 
performed using the Navy version of AQAM and 1975 activity at 
NAS Miramar as a representative data base (Ref. 9). 

The purposes of the present study were (1 ) to update the 
data in the Source Inventory program of AQAM in order to rep- 
resent 1978/1979 operations at NAS Miramar and (2) to compare 
the predicted and measured levels of pollutant concentrations 
for the purpose of validating the Short Term program of AQAM. 

A necessary component of the validation effort was the con- 
ducting of an updated model prediction sensitivity study. 

II. OVERALL MODEL VALIDATION EFFORT 

The Navy version of the AQAM model validation effort was 
initiated by the Naval Air Propulsion Center (NAFC). NAPC 
provided the funding and necessary program coordination as 
well as technical assistance in selection of the monitoring 
site locations and the required data acquisition. NAS Miramar 
was chosen because it had the largest number of flight opera- 
tions of any NAS and because it had been used in previous work 
performed by the Naval Postgraduate School in developing the 
Navy version of AQAM. 

The overall objectives of the NAPC program were to: 

a. validate the AQAM model, 

b. document the effects of aircraft operations on 
air quality, and 
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c. assess the possibility of using AQAM (as an alter- 
native to an expensive monitoring program) to 
determine the effects of aircraft operations on 
air quality at other NASs (Ref. 11). 

The program was divided into two related parts. The first part 
was to consist of a one year continuous monitoring study. In 
late 1979 - early 1980 air quality was being measured 24 hours 
a day using an automated data acquisition system. This effort 
was directed primarily at objective (b) noted above. The second 
part was to consist of two special studies, each one-week in 
duration. The latter studies were to be intensive in nature 
with detailed operational, meteorological and pollution concen- 
tration data being collected. These studies were to be directed 
primarily at objectives (a) and (c) above. The first special 
study took place in August 1979 and data received from that week 
were used in the model validation discussed herein. The second 
special study was scheduled for the spring of 1980 but was sub- 
sequently cancelled. The two periods were chosen to occur during 
distinctly different meteorological conditions, especially lid 
height and stability category. Organisations involved in the 
special study and individual responsibilities of each included: 

a. Northrup Services Incorporated (NSI) contracted 

by EPA: Air quality monitoring and data reduc- 

tion to provide hourly averaged pollutant levels. 

b. Pacific Missile Test Center (PMTC) : Meteorological 

measurements and data reduction to provide hourly 
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averaged weather conditions throughout the 
receptor grid. 

c. NAPC/NPS: Aircraft activity monitoring. 

d. NPS: Reduction of aircraft activity data for in- 

put into A QAM, model predictions using items b. 
and c. above, and comparison of predictions with 
measured values. 

III. NAS MIRAMAR INTENSIVE DATA ACQUISITION 

Planning the special study for validation of AQAM began 
■with identifying both the emittants to be monitored to best 
characterize dispersion and, as previously mentioned, locating 
appropriate monitoring stations. 

The major pollutants in aiccraft engine exhausts include 
particulates/smoke (?T), carbon monoxide (CO), unburned hydro- 
carbons (HC) and oxides of nitrogen (NOX). The relative amounts 
emitted depend primarily upon the engine thrust setting. In 
addition, sulfur oxide emissions (SOX) are often significant 
from industrial and domestic furnaces. Therefore, CO, HC, NOX, 
PT and SOX were selected as the pollutants to characterize 
emissions of both aircraft and airbase related activity. Figure 
1 identifies the grid system used to locate the receptors in * 
AQAM. The grid spacing was 1 km. and the x-y coordinates 
varied from (0,0) to (21+, 15 ) representing 1+00 separate recep- 



tor locations 
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GRID IDENTIFICATION AT NAS MIRAMAR 



Locating measuring stations where continuous-air-monitor- 
ing instruments would be placed was of prime importance in the 
validation effort. The behavior of the model predictions at 
a particular receptor will depend to a great extent on its 
location relative to numerous sources throughout the receptor 
grid, especially those located upwind. To validate the model, 
it was important to compare air quality samples at locations 
where the airbase and aircraft contributions were large rela- 
tive to background levels of pollution. Ultimate placement of 
the stations assumed a dominant wind from the WNW (292°) as 
advised by PMTC. 

Up to 12 special receptor locations can be input to the 
Short Term program. Special receptor locations were assigned 
to each of the four pollution monitoring stations as indicated 
in Table I. They are also identified in Figure 1 . 



TABLE I 

MONITORING STATION LOCATIONS 



TRAILER NUMBER 


GRID COORDINATES 


SPECIAL RECEPTOR 
NUMBER IN A QAM 


1 


1 0 . 01 , 8.24 


401 


2 


10.52, 8.46 


402 


3 


11 . 24 , 8.35 


406 


4 


12 . 82 , 7.31 


410 



The intended use of trailer 1 was to determine background 
levels of pollution upwind of aircraft/airbase sources. 

Trailer 2 was located just downwind of the hot refueling site. 
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Trailer 3 was situated just upwind of the hot refueling pits. 

It was also downwind of the hot refueling area. Trailer i; 
was located well downwind at the outer boundary of MS Miramar. 

During the planning stage, NSI made equipment preparations 
for each trailer site for the air quality monitoring experi- 
ment. PMTC analyzed the meteorological history for the San 
Diego area to determine the best time period for the special 
study. Optimum weather conditions for validation were con- 
sidered to consist of a moderate wind coming from the 290 
degree direction, a Turner stability category of 2-3, and a 
moderate lid height (mixing layer depth) of lj.00-500 meters. 

It was desirable to have relatively constant weather condi- 
tions for the week of intensive data acquisition. This would 
allow the dispersion model to be validated with multiple tests 
in which aircraft operations varied but weather remained approx- 
imately fixed. The week of 1 -? August 1979 was chosen as the 
most feasible for meeting these objectives for the first in- 
tensive study. 

Operating procedures for the week proceeded on a previously 
planned routine. Specific tasks performed by NSI (pollution 
monitoring) and PMTC (meteorological monitoring) will be pre- 
sented by those activities under separate cover. NFS and NAPC 
personnel monitored the detailed aircraft activity in accor- 
dance with the time schedule listed in Table II. 
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TABLE II 



AIRCRAFT ACTIVITY MONITORING TIMES (LOCAL) 



1 AUG 



1300 - 1600 



2 AUG 



1000 - 1230 
1 i+oo - 1700 



3 AUG 



0800 - 1230 



6 AUG 



0900 - 1230 
1 330 - 1 630 



7 AUG 



0830 - 1100 

1 330 - 1 630 



Observation of aircraft activity was performed/recorded 
from three locations — the control tower, the hot refueling 
site (octagon) and the refueling pits. 

The functions performed in the control tower involved 
(1 ) timing the sequences of every aircraft on departure from 
initial startup to takeoff, (2) timing the sequences of every 
aircraft on recovery from entry into the airport traffic area 
(defined here as having a three-mile radius) to landing and 
taxi to the refueling area. Also, the parking areas and taxi- 
ways used by each aircraft and the type of landing performed 
(VFR, IFR) were monitored. Data sheets used to record the 
aircraft activities observed from the control tower are pre- 
sented in figures 2 and 3* 

Data collected at the hot refueling sites (octagon) and 
refueling pits included time-in-mode, amount of fuel taken, 
and aircraft type (see data sheets in figures if. and 5)« 
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TAKEOFF DATA SHEET 



DUTY RUNWAY WIND 

register/time 

Side number Aircraft type 



Parking area 


Commence sequence 


0/0 


Start complete 


/ 


Taxi complete 


/ 


(holding at runway 
(engine check comp 


7 " 

lete) 


Takeoff complete 


/ 



EVOLUTION (check one) 

Takeoff and deparc area 

FCLP number 

Touch and go number 



FIGURE 2 



LANDING DATA SH5ST 
(full stop landings only) 

DUTY RUNWAY WIND TIME 



register/time 

Side number Aircraft type 

Commence sequence 0/0 

(enter break or 3 mi. on IFR approach) 

Landing complete^ 

(clear of runway ) 

Taxi complete 

(pits/hot refuel holding area) 

Fuel commence 

(enter pits/hot refuel area) 

Fuel complete 

(depart pits/hot refuel area) 

Shutdown 

Parking area 

(hot refuel aircraft only) 



FIGURE 3 
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TIKE 



HOT REFUEL SEQUENCE DATA SHEET 
( OCTAGON ) 



A/C Type (circle one) 



Side number_ 

Arrival time 
at 

holding area 

Arrival time 
into 

octagon 

Departure time 
from 

octagon 

Pounds fuel 

received 

Fuel spilled yes/no 



P-5 

P-8 
P-1 ix 

A«V 

3-2 

Transient 



FIGURE ij. 
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PIT REFUEL SEQUENCE DATA SHEET 



TIME 



A/C Type (circle one) 



Side number_ 

Arrival time 
at 

holding area 

Arrival time 
into 

refuel pit 



P-4 

P-5 

P-8 

P-14 

A-4 

E-2 

Transient 



Shutdown 

(circle one) 

Pounds fuel 

received 

Fuel spilled ?es/no 



Hot refuel 



Departure time 
from 

refuel pit 

Pounds fuel 

received 

Fuel spilled yes/no 



FIGURE 5 
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The aircraf t/airbase operational data that were collected 
were used as input to the Source Inventory program. Air 
quality measurements (by NSI) and meteorological data (by 
FMTC) were also being collected during the entire period of 
observation. 

IV. A QAM MODIFICATIONS AND SENSITIVITY STUDY 
A. MODEL MODIFICATIONS 

In order to perform a model validation, the data input to 
the Source Inventory program must reflect, as closely as pos- 
sible, conditions and emittant sources as they exist at the 
time of validation. Therefore, one of the purposes of this 
study was to update the data in AQAM to represent 1973/1979 
operations at NAS Miramar. 

Changes made to the input routines of the AQAM program 
included data input on the S-2 aircraft — an addition at NAS 
Miramar since 1975. Parking area coordinates, taxiway usage 
and aircraft landing and take-off operational cycle time-in- 
mode (LTO) were all modified to accept S-2 aircraft activity. 
All data were input in accordance with guidelines stipulated 
in Refs. 7-9 and 12. Averaged meteorological data were 
changed to reflect 1973 figures. The annual amount of air- 
craft activity for 1973, including arrivals, departures, touch- 
and-go cycles, and FCLP’s was entered according to aircraft 
type. The specific parking areas and taxiways used by each 
aircraft were modified. Other emissions information 
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(specifically; fuel spillage, training fires, environ land 
use area factors, and off base vehicle miles per year) was 
either added or updated. Airbase, non-aircraft activity 
modifications included changes in test cell and run-up stand 
usage. 

3. SENSITIVITY STUDY PARAMETERS AND PREDICTIONS 

With the update completed and reflecting conditions as 
they existed at the time of the first intensive study, an 
investigation was performed to determine the sensitivity of 
the model predictions to meteorological and operational con- 
ditions anticipated for 1-7 August 1979 (special study). Sen- 
sitivity results indicate under what conditions and at what 
receptor locations the model can best be validated. In addi- 
tion, these results are needed before conclusions can be drawn 
from the comparison of measured and predicted pollution levels. 
In a model validation effort, predicted concentrations are 
compared to measured values at specific receptor locations. 

When making these comparisons it is necessary to know how sen- 
sitive the model predictions are to the uncertainties in the 
specified meteorological and operational input data. For 
example, stability category is normally specified as an inte- 
ger value between one and six; if the hourly averaged value 
can only be specified as two or three, what effect would this 
variation have on the model predictions? In addition, it is 
necessary to know whether the monitoring stations are located 
in regions where there are large horizontal gradients in pol- 
lution concentrations. 
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Twelve special receptors were used to examine the sensi- 
tivity of predicted pollution levels in the vicinity of the 
four monitoring stations to various meteorological conditions 
and model parameters. A previous model sensitivity study had 
been conducted by Netzer (Ref. 9) using 1975 operational data 
and different nominal meteorological conditions. Table III 
describes the special receptor locations used in A QAM for both 
the sensitivity study and the validation effort. Locations 
relative to runways, taxiways and refueling areas are depicted 
in Figure 1 . 



TABLE III 

SPECIAL RECEPTOR LOCATIONS 

A QAM RECEPTOR NUMBER DESCRIFTION/LOCATIOH 



401 


trailer 1 




402 


trailer 2 




k03 


100 m downwind 
trailer 2 


of 


404 


100 m crosswind 
of trailer 2 


( south) 


k05 


100 m crosswind (south- 
east) of trailer 2 


406 


trailer 3 




k07 


100 m downwind < 
trailer 3 


of 


403 


100 m crosswind 
of trailer 3 


( south) 


409 


approach end of 


runway 1 


410 


trailer 4 




41 1 


500 m upwind of 

4 


trailer 


CM 
* — 


100 m crosswind 
of trailer 4 


(north) 
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In order to perform the sensitivity study it was neces- 
sary to establish a reference or nominal case meteorologically 
and operationally* The anticipated weather conditions for the 
intensive study period, listed in Section III, were used as 
the reference weather. Meteorological parameters were varied 
independently, with aircraft activity kept constant. Table 
IV indicates the meteorology data input for each of nine com- 
puter runs. 



TABLE IV 

METEOROLOGY FOR SENSITIVITY STUDY 

Rim Turner Wind Wind Temperature Lid 

Number Stability Speed Direction /O n \ Height 

(m/s) (deg) { ^ (m) 



1 2 i+. 1 2 290 80 I4.OO 
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Run number 1 was the reference case. The ambient air temp- 
erature was not varied because previous results (Ref. 9) had 
shown it to have little effect on predicted pollution levels. 

The aircraft activity data input to the Source Inventory 
program were representative of one hour of daytime flight 
operations. In addition, airbase and environ sources were 
kept constant with updated 1979 data. In the normal mode of 
utilization of AQAM, annual totals are input and frequency 
factors are used to determine the total operations in any one 
month, week, day, and hour. For the present effort, the 
"desired 11 one hour input data had to be scaled up to annual 
operations in order that the Short Term and Source Inventory 
programs would function properly. The "scale-up" factor used 
was : 

12 hr/aay x 31 day /mo (Aug) x 12 mo/yr = kii6k hr/yr (1 ) 

(12 hr/day represents no night operations) 

Table V presents the aircraft activity values which were held 
constant for the entire sensitivity study. 
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TABLE V 



AIRCRAFT ACTIVITY FOR SENSITIVITY STUDY 



\ riOUn OrhitlAl J_Ur< o 


AIRCRAFT 


ARRIVALS 


DEPARTURES 


TOUCH 

& 

GO'S 


VFR 

ARRIVALS 


FCLP'S 


F-4 


3 


3 


2 


2 


6 


F-3 


1 


1 


1 


1 


0 


E-2 


1 


1 


1 


1 


0 


p-i u. 


3 


3 


2 


2 


6 


A-U. 


2 


2 


1 


1 


0 


P-5 


1 


1 


0 


1 


0 


TRANSIENT 


1 


1 


0 


0 


0 


H-3 


0 


0 


0 


0 


0 



1 YEAR OPERATIONS 



AIRCRAFT 


ARRIVALS 


DEPARTURES 


TOUCH 

2c 

C-O's 


VFR 

ARRIVALS 


FCLP'S 


P-4 


13392 


1 3392 


8928 


8925 


26 784 


F-3 


144.61; 


4464 


4464 


4464 


0 


E-2 


414.614. 


4464 


4464 


4464 


0 


F-14 


13392 


1 3392 


8928 


3928 


26784 


A-4 


3928 


8928 


4464 


4464 


0 


P-5 


4464 


4464 


0 


4464 


0 


TRANSIENT 


44.61; 


4464 


0 


0 


0 


H-3 


0 


0 


0 


0 


0 



As explained in Section I, the results from the Source In- 
ventory program are used along with the meteorological data as 
input to the Short Term program. Output from the Short Term 
program was arranged in seven tables. Four tables consisted 
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of pollutant levels in micrcgrams per cubic meter from environ, 
airbase, aircraft and total sources at all specified grid re- 
ceptors. Each table listed, for all receptors, the receptor 
number and x-y coordinate location, and the concentrations for 
all five pollutants. The remaining three tables expressed the 
same results in terms of fractions of the total emissions from 
environ, airbase and aircraft sources. 

The receptors of interest in the sensitivity study were 
the twelve special receptors (1±01 -i|.1 2 ) and that one where the 
maximum concentrations existed. 

To compare the expected effects of the meteorological 
variables on the predicted ground level (z=0) concentrations, 
the Gaussian dispersion formula for point sources can be used 
(Ref. 13). 



:(U,y,z=o jh)» 

where : 



7rcr y TTS 



exo 



2 V 



exo 



_L_ ( JL)2 

2 v ' 



( 2 ) 



v* = 



concentration. 



g/n 






Q = uniform emission rate, g/sec 



<r 



y* 



standard deviations of plume concentrations 
in the horizontal and vertical directions 
respectively, m 



U = mean wind speed, m/sec 



H = initial plume height, m 
y = 0 along plume centerline 



When vertical diffusion is limited by a stable layer at 
height h, . , the diffusion equation is modified as follows: 
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(3) 



: a 

vTu <r h, .,17 

y 



exp 




y 



?or infinite line sources Turner (Ref. 13) utilized: 



X(x,y.,z=0;H) = 



lQ 



sin $ fZv* U 
z 



exp 



_ 1 _ 

2 



^> 2 



( 1 +) 



■where : 

q = source strength per unit distance, g/sec-m 
= angle between line source and wind direction, 

k5°< < 90° 

Major variations of the Short Term program predictions under 
different meteorological conditions should follow equations 
(2), (3 ), or (I4.), depending upon the receptor location rela- 
tive to the dominant emission sources (Ref. 9), 



C. EFFECT OP METEOROLOGICAL PARAMETERS ON MAXIMUM RECEPTOR 

CONCENTRATIONS 

Table VI presents the predicted maximum concentrations 
of four of the five pollutants and the location of each for 
the reference case. Also shown are the maximum predicted CO 
and PT from aircraft sources for each of the other conditions 
investigated. The meteorological variable is listed in each 
case . 

The reference case indicated that the maximum contribu- 
tions from the environ sources occurred south of the airbase 
(at receptors (9,2) and (11,2)). However, high levels of environ 
pollution (background) also were predicted to occur throughout 
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MAX J HUM aHOmPTOR CONCENTRATIONS 
(CHID LOCATION, CONCENTRATION (/>c;in/m^ ) , FRACTION OF TOTAL 
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the airbase 



On the airbase, the contribution from airbase 



sources was generally negligible, whereas aircraft sources of 
FT were dominant. Maximum concentrations from aircraft sources 
occurred for CO and PT at receptor (11,3), near the intersection 
of the runways. This was generally the case for all the con- 
ditions investigated. 

1 . Stability Category 

Increasing the stability category (more stable condi- 
tions) decreases <r and ar , and therefore should increase the 

y z 

predicted ground level concentration along the wind vector 
downwind of the source (see equation (2)). At the peak concen- 
tration receptors (Table VI), which are necessarily near the 
plume centerline, the increase in stability category increased 
the concentration and shifted the maximum concentration receptor 
downwind. 



2. Lid Height 

As a plume develops downwind of a source it will spread 
in a vertical, as well as horizontal, direction. The ground 
and lid height (elevated inversion layer) act as reflectors of 
the plume. Increasing the lid height would decrease the con- 
centration only at receptors which are far enough downwind from 
the source for reflections to occur (see equation (3)). For 
the maximum receptor location (11,3), lid height had negligible 
effect on the predicted aircraft concentrations of CO and PT (Table VI) 
since it was located near the major aircraft sources. However the aircraft 
contributions increased with lid height since concentrations from environs 



decreased . 
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3. Wind Speed 



Increasing the wind speed should decrease predicted 
concentrations along the plume centerline for a single source 
(equations (2), (3) and (1|)). Thi3 behavior was apparent for 
the maximum concentration receptors (Table VI, run nos. 6, 1 
and 7 ) . 

!(.. Wind Direction 

Changing wind direction changes the orientation of the 
plume dispersion. As a result, the receptor where concentra- 
tions were a maximum from aircraft sources was predicted to 
shift to receptor (10,8) when the wind direction became 270° 
(Table VI, run no. 8). 

D. EFFECT OF METEOROLOGICAL PARAMETERS ON CONCENTRATIONS 

AT SPECIAL RECEPTORS 

Short Term output for each of the nine sensitivity runs is 
presented in Appendix A for the special receptors. The refer- 
ence case (run no. 1 ) output includes receptor concentrations 
for environ, airbase, aircraft, and total sources in //gm/m^ as 
well as fractional values for aircraft sources. Receptor con- 
centrations for aircraft sources (run nos. 2-9) are included 

•3 

in and fraction of total. In order to visualize varia- 

tions in pollutant concentration, che overall grid system was 
mapped with contour levels for the sensitivity study in Appen- 
dix 3. Contours for the reference case are included for CC 
and PT concentrations from airbase, aircraft, and total sources. 
Contours for run nos. 2-9 are included for CO and PT concentra- 
tions from aircraft sources. 
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Tables VTIa-d summarize the special receptor concentra- 
tions of CO and PT for each of the nine sensitivity runs. In 
general, the comments relating to the maximum receptor concen- 
trations pertain to the special receptor concentrations. From 
a modeling standpoint special receptor lj.01 (trailer 1 ) proved 
to be well located for the purpose of measuring background 
pollutants. As can be seen in Tables Vlla-d, very little CO 
and PT due to aircraft exist at receptor lj.01 . When finite 
values did occur (run nos. 2, 3> 7 > 9 and 9) they resulted 
from the aforementioned method of projecting area sources (in 
this case — the hot refueling area) upwind to pseudo -point 
sources. 

1 . Stability Category 

An increase in stability category increases the down- 
wind concentration along the plume centerline from a single 
source since the plume spreads more slowly. Table Vila indi- 
cates that the area around trailer 1 (receptors I4.02-I4.05) 
receives emittants from multiple sources since the concentra- 
tions of CO and PT first decreased and then increased with in- 
creasing stability category. These receptors are also located 
very near large sources. 

CO and PT concentrations around trailer 3 (receptors I4.O6- 
I4.O8) were significantly higher than those around trailer 2 due 
to the effect of an increased number of plumes overlapping 
downwind. Some multiple/near source effects were also evident 
at this location. The receptor concentrations around trailer 
I4. (receptors l4.IO-lj.i2) changed only slightly with variations 
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TABLE Vila 
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TABLE VI lb 
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in meteorological conditions due to the large downwind distance 
from the primary sources. Concentrations at receptor lj.09 were 
high as expected due to its close proximity to runway and taxi- 
way line sources. 

2. Lid Height 

At trailers 2 and 3 lid height had no effect (Table 
Vllb). This was expected since these locations are very near 
the sources of pollution. At trailer ij., which is far down- 
wind, increasing lid height decreased concentrations. 

3. Wind Speed 

As indicated in Table Vile, an increase in wind speed 
decreased the concentration downwind at trailers 3 and Ij.. 

Again, however, at trailer 2 the behavior was more random. 

Ij.. Wind Direction 

Changing the wind direction from the reference 290° 
to 31 0° (run no. 9) resulted in the expected reduction in air- 
craft CO and ?T at trailers 2 and 3 (Table Vlld). In this case 
the plumes from the major upwind aircraft sources miss receptors 
1+02 and 14.06 . However, when the wind direction was changed to 
270° (run no. 8), the concentrations increased significantly. 
This indicates that trailer 2 was apparently outside the plume 
from the hot refueling area when the wind was from 290°. 

Further evidence of this was that receptors iiOl; and Il 05 (cross- 
wind to I4 02) had significantly higher concentrations than re- 
ceptors ij.02 and I4O3. 

The trailer k receptor exhibited an increase in concentra- 
tion with an increase in wind direction. This was expected 
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since most aircraft source plumes (including the maximum re- 
ceptor location at coordinate (11,8) are located upwind of 
trailer l)., from the 290°-310^ direction. 

5 . Special Receptor Locations 

As discussed above, for model validation efforts it is 
necessary to know whether the monitoring stations are located 
in regions where there are large horizontal gradients in pol- 
lution concentration or where the concentrations are very 
sensitive to the specified hourly-averaged meteorological con- 
ditions. Table VIII presents a summary of the effects of 
distance from the monitoring stations on the predicted pollu- 
tion concentrations. Concentrations are presented for each of 
the nine cases for conditions 100m downwind and 1 OCm crosswind. 

As a receptor is moved toward a specific plume centerline, the 
concentration would increase. When a receptor is located down- 
wind from several sources, horizontal movement of the receptor 
may increase or decrease the receptor pollution level, depending 
on the multiple plume effects. 

Increases in concentration varied by factors of two to 
sixteen at trailers 2 and 3 for the reference case as a result 
of moving the receptor 100m downwind or closer to plume center- 
line. No appreciable horizontal gradients in concentration 
existed around trailer ij.. In almost every case (variation of 
meteorological parameters), concentrations increased as expected, 
since the receptors were moved closer to the centerlines of the 
major aircraft-related plumes for the 290° wind. In run no. 8, 
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TABLE VIII 

DIFFERENCE FACTORS IN SPECIAL RECEPTOR CONCENTRATIONS 
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where the wind direction was changed from 290° to 270°, the 
concentration at the 100m downwind location decreased at 
trailer 2. 

These results again indicate that comparison between 
measurements and predictions will be most difficult at trailer 
2. Not only do multiple plume effects and the close proximity 
to ground aircraft sources cause unusual variations in con- 
centration with changing meteorology but also the horizontal 
gradients are quite large. 

E. EFFECT OF SPECIFIED AREA SOURCE SIZE ON RECEPTOR 

CONCENTRATIONS 

’//hen large sources are input into AQAM they are normally 
modeled as area sources. The dimensions of the area sources 
have to be specified and some judgement is required to pick 
the most representative dimensions of these "uniform concen- 
tration sources." To determine what effect the specified 
size of aircraft area sources had on concentrations at various 
receptors, the lengths of the sides of three prime sources 
were both increased and decreased by forty percent. The 
specific sources included the hot refueling area, the hot re- 
fueling delay area and the pit refueling delay area. The length 
of the sides of each area source in the reference case was 500 
meters. This length was changed to 300 meters and then to 700 
meters . 

Increasing the size of an area source effectively moves 
the pseudo-upwind point source further upwind. Keeping the 
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emittant3 and meteorology constant, the plume would spread 
at the same rate. At a specific receptor, the concentration 
can increase or decrease, depending on its location relative 
to the area sources. For this study, the variations in con- 
centrations at trailers 2, 3> and if. never exceeded six percent. 

?. VARIATION 0? JET PENETRATION LENGTH AND HORIZONTAL AND 

VERTICAL DISPERSION PARAMETERS 

In AQAM, turbojet exhausts during taxi and takeoff are 
treated as finite line sources. Initial line source dimen- 
sions and locations have to be specified and these are some- 
what arbitrary. Currently in AQAM the jet is assumed to 
''penetrate IliO meters' 1 (i.e., approximately 1 2x0 jet diameters) 
before coming to rest relative to the ambient air. Default 
values for the line source cross-sectional size are 3m high 
by 20m wide. No plume rise i3 considered to occur. These 
line sources are then treated as pseudo-upwind lines which 
disperse in a Gaussian manner with the same empirical disper- 
sion parameters as used for elevated point sources. 

In a recent study at the Naval Postgraduate School (Ref. 

U|.) jet characteristics were measured in a simulated, neutrally 
stable atmosphere. It was found that jet penetration length 
was considerably less than llf.0 jet diameters; being more nearly 
35 jet diameters. Initial plume dimensions were found to vary 
significantly with jet orientation to the ambient wind direc- 
tion and some plume rise was observed. Jet dispersion rates 
were also found to be more rapid than currently used in AQAM. 



In order to determine whether the above findings would 
have any significant effects on the predicted concentrations 
from aircraft sources, AQAM was modified in sequential steps 
as follows: 



(1) decrease the jet penetration length from 140 
to 35 meters. 

(2) step (1) and specification of initial aircraft 
line source (taxiway and runway) dimensions as 
a function of orientation to the wind (per fig. 
40, Ref. 14). 

(3) steps (1) and (2) and decrease the stability 
category by one to increase the jet plume 
spreading rate. 

Decreasing the penetration length was found to have little 
effect. This was somewhat expected since the aircraft line 
sources at NAS Miramar have lengths up to 3.7 km. The reduc- 
tion in jet penetration length was but three percent of the 
longest line source. In step (2) the angle of incidence 
formed by the wind with each line source was determined, and 
using the a^ r and a „ versus angle of incidence relationship 
determined by Brendmoen and Netzer, new horizontal and vertical 
dispersion parameters were input to the Short Term program. 

In general, the changes involved increases in initial line 
source dimensions. At the maximum concentration receptor and 
at trailers 3 and 4, a nominal reduction in concentrations of 
up to a maximum of 16 percent was predicted. 
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In step (3) the above changes were kept in A QAM and the 
stability category was decreased from 3 to 2 (more unstable 
conditions). Output indicated a decrease in concentration of 
up to a factor of two at the maximum concentration receptor 
and at trailers 3 and Ij.. It should be noted that in its pre- 
sent form A QAM only allows variation of stability category for 
all dispersions as opposed to variation of aircraft sources 
alone. This decrease was expected as previously determined in 
the meteorological sensitivity study. 

G. CONCLUSIONS 

Stability category and wind speed were the two meteoro- 
logical parameters that most affected maximum receptor concen- 
trations. Model predictions will therefore be most sensitive 
to uncertainties in the hourly-averaged values of these para- 
meters which are input into AQAM. Wind direction had a large 
effect on the concentrations at trailer 2. Trailer 2 is 
apparently located in an area where large horizontal gradients 
of pollutant concentrations exist, i.e., near the edges of the 
plumes from large aircraft sources. 

Trailer 1 appears to be a good location for measurement of 
background pollution levels. 

Variations in aircraft area source sizes did not appreciably 
affect concentration levels at specific receptors. 

Variations of the specified jet penetration length and 
initial horizontal and vertical dispersion parameters of air- 
craft exhaust plumes during taxi, takeoff and landing modes 
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changed concentrations by a maximum of only 16 percent. The 
data of 5rendmoen and Netzer (Ref. lij.) indicated that turbo- 
jet exhausts spread more rapidly than elevated point sources. 
This result, when incorporated into A QAM, significantly 
affected predicted concentration levels (by a factor of 2) at 
the monitoring trailer locations. 

% 

V. COMPARISON OF A QAM PREDICTIONS 

WITH 

DATA FROM THE INTENSIVE STUDY 
A. VALIDATION REQUIREMENTS 

As previously stated, model validation consists of compar- 
ing predicted hourly- averaged pollutant concentrations to 
hourly-averaged measured values at specific receptor locations. 

A determination of model accuracy must be made within the con- 
text of the accuracy of the input operational data and of the 
hourly-averaged meteorology and measured concentrations. It 
is important to note that although the meteorology and pollu- 
tant concentrations may be constantly varying, only hourly- 
averaged values are used. Comparisons between measured and 
predicted concentration values in areas where large horizontal 
gradients exist (trailer 2) are likely to exhibit widely-varying 
results. Because of these factors, a need exists for a vast 
amount of accurate data with which to conduct model validation. 

Prior to the comparison of measured and predicted values, 
background levels/local air quality must be determined in order 
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to be able to separate the contributions of 1 aircraft, airbase 
and environ sources throughout the receptor grid. The Source 
Inventory program allows for input of environ sources. If 
these data are not available, approximate inputs can be in- 
cluded through the use of land-use factors. The factors (Ref. 
12) distinguish between city center, urban, rural, park areas, 
etc. Input for off-base line sources (roadways) requires 
appropriate vehicle mileage and speed values. The selection 
of appropriate land-use factors used in this study was some- 
what judgemental. The roadway line source values used were 
based on actual average daily traffic volumes for 1973 as pro- 
vided by the Comprehensive Planning Organization of the San 
Diego Region. One method for determining actual concentra- 
tions from aircraft/airbase sources is to subtract values from 
an upwind measurement (i.e., trailer 1 data) from values 
obtained at each of the other special receptors. 

Comparison of weekend measured data at each special recep- 
tor with weekday data should also provide a good indication of 
background/environ pollutant levels due to the reduction in 
aircraft activity at NAS Miramar on weekends. The measured 
data indicated a wind speed varying from calm to five knots 
on Saturday and Sunday approximately 90^5 of the time. The 
wind direction also varied up to 1 80° throughout the two-day 
period. This slight-to-stagnant air motion apparently caused 
an accumulation of pollutants at NAS Miramar from environ 
(local San Diego) sources. Unfortunately, this behavior 
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invalidated any comparison between weekday and weekend concen- 
trations for the purposes of validating weekday background levels 
on the airbase. Therefore, a need exists for additional weekend 
data when the meteorological conditions are more representative 
of those experienced during the period of intensive measurement. 

B. DATA REDUCTION AND MODEL INPUTS 

Measured data for CO, NOX and THC were provided by NSI in 
parts per million (ppm). Comparison of these values to AQAM 
predictions requires conversion to micrograms per cubic meter 
(^gm/m^) . An accurate conversion exists for CO under standard 
conditions; 1111.11 x ppm CO = > mgm/m^ CO. The most often used 
conversion for NOX is based upon N0 9 : 2000 x ppm NOX =yt 4 gm/m^ 

NOX. Measured data were obtained for THC and CH^ . CH^ often is 
the major portion of THC concentration in urban atmospheres of 
North American latitudes. Typical concentrations are 1.25-1.5 
ppm (Ref. 6). The CH^ conversion is 666.67 x ppm CH^ = .^gm/m^ 

CH^ . However, aircraft generated hydrocarbons may be significant- 
ly heavier than CH^. The only PT data available were measured 
by a nephelometer in terms of the scattering coefficient, b (bscat) 
Air samples were also taken to determine total particulates (TP) , 
but the data were invalidated as a result of a filter preparation 
error by contractors at U. C. Davis. This loss of the TP data 
severely affected the model validation effort. Particulate con- 
centrations on the airbase had the best chance to be dominated 
by aircraft sources and therefore provided one of the best means 
for comparing predictions with measurements. For the bscat data, 
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an average conversion factor was employed (Ref. 15); 46.15 x 
Neph (bscat) ~ ^gm/m 0 PT. The latter conversion factor may be 
in considerable error. 

The AQAM model was run over ten one-hour time periods as 
listed in Table IX. The type of aircraft activity varied con- 
siderably throughout the ten AQAM runs. When different from 
normal operations, remarks of the activity are included in Table 
IX. The chosen periods of time were primarily in the afternoon 
when the wind speed and lid height are greatest. 

Figures 6a-e present the meteorological conditions at NAS 
Miramar (obtained from NAS, Pt . Mugu investigators) for the days 
of intensive measurements and detailed observation of aircraft 
activity. The values are hourly- averaged and plotted over the 
1000-2000 time period for each day. All weather conditions were 
averaged over the applicable time periods shown cross-hatched 
in Figures 6a-e. 

Runs 5 and 9 were performed to determine ivhether or not 
transit time of emittants affected predicted concentration levels 
relative to runs 4 and 8. A fifteen minute emittant travel time 
was chosen due to the wind speed and average distance from source 
to monitoring station. It should be noted that the final runs 
(6-7 August) had significantly higher wind speed, temperature and 
lid height. This variation in meteorology was not anticipated 
and was somewhat undesirable from a model validation viewpoint. 

Due to variations in meteorology within the calculated 
dispersion times, it is generally agreed that the values of 
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METEOROLOGICAL CONDITIONS AT NAS MIRAMAR VS. TIME OP DAY (2 AUG) 
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METEOROLOGICAL CONDITIONS AT NAS MIRAMAR VS. TIME OF DAY (3 AUG) 
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METEOROLOGICAL CONDITIONS AT NAS MIRAMAH VS. TIME OP DAY (6 AUG) 
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METEOROLOGICAL CONDITIONS AT NAS MIHAMAR VS. TIME OF DAY (7 AUG) 



CT and d_ cannot be more accurate than a factor of 2. In 

y 

addition to this uncertainty, model predictions are sensitive 
to the average meteorology used as input as discussed above. 

For example, consider the data presented in Fig. 6a for the 
period 1300-1400 hrs . During this period the wind direction 
changed from 270° to 300° and the stability category changed 
from 1 to 2. The lid height and wind speed were steady. Values 
employed for wind direction and stability category for this period 
(Table IX, run no. 1) were 290° and 2, respectively. The sensi- 
tivity study of section IV has shown that a decrease in wind direc 
tion of 20° and a decrease in stability category from 2 to 1 can 
increase the predicted concentrations at trailer 3 by factors of 
1.5 and 1.3 respectively. Thus, measured data and predictions 
could be different by a factor of approximately 2 due to uncertain 
ty in model meteorological input alone. 

C. DISCUSSION OF RESULTS 

This discussion is divided into four sections -- one for 
each of the pollutants measured. The included figures are scatter 
plots of measured CO, NOX, THC and nephelometer readings versus 
predicted concentrations. The diagonal lines drawn in these fig- 
ures enclose predictions that are within a factor of two of the 
corresponding measurements. These lines were found to enclose 
greater than 50% of all the plotted points. Much of the measured 
data were invalidated by NSI and were therefore not available for 
plotting. This is the reason for the differences in numbers of 
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plotted points from graph to graph and was a major limitation 
in the validation effort. In fact, only approximately 451 of 
the pollution data taken during the intensive study were con- 
sidered acceptable. None of the data taken during the first 
three days of the intensive study met the validation criteria 
in the EPA quality plan. However, those data had to he used in 
order to have even a minimum of data to compare with model pre- 
dictions . 

Variations in predicted pollutant concentrations over the 
airbase were mapped with contour levels for the intensive study 
and are presented in Appendix C. Contours for run no. 4 (2 Aug, 
1500-1600) are included for CO and PT concentrations from air- 
base, aircraft and total sources. Contours for the other nine 
runs are included only for CO and PT concentrations from air- 
craft sources. 

1 . Carbon Monoxide (CO) Emissions 

A comparison of the CO emitted each weekday with the 
CO emitted during the weekend (period of reduced aircraft 
activity) was performed to better determine the CO background 
level. It was found that on Saturday afternoon the level was 
higher than that on Monday by a factor of two, possibly due 
to heavy traffic conditions on the surrounding roadways. 

Also, on Sunday, when the winds were mostly calm or from the 
south, a high level of CO was measured at trailer 1. As pre- 
viously stated, weather conditions for the weekend during the 
period of intensive measurement were not representative of 
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weather conditions during the weekdays. Therefore, no con- 
clusions could be drawn from this comparison regarding the 
validity of using trailer 1 measurements as indicators of 
background CO levels. 

Figure 7 indicates that measured concentrations agreed 
with predicted total concentrations within a factor of two at 
trailers 1 and 4. The agreement was within a factor of approx- 
imately three for trailer 3 data. (No measured CO data for 
trailer 2 was available during the ten one-hour time periods 
used in this study.) However, the good agreement may be chance 
since the environ input (land-use factors, vehicle mileage data, 
etc.) was only estimated. In other words, what if the high 
levels of CO concentration at trailer 1 were due to aircraft, 
but the model did not have properly input aircraft operations 
or did not correctly determine dispersion rates? AQAM predicted 
that the CO concentration due to aircraft at trailer 1 was 
essentially zero. To check this, the Source Inventory program 
was modified so that all aircraft climb angles on takeoff were 
decreased. This maximized the near ground emissions from air- 
craft in the area near trailer 1. This change had no effect on 
CO at trailer 1. Also, the sensitivity study discussed above 
indicated no effect from increasing the hot refueling area and 
hot refueling delay area source sizes. In other words, some 
inaccuracy in aircraft source specification near trailer 1 would 
not cause increased concentrations at that receptor. Therefore, 
it appears to be a valid assumption that trailer 1 was a good 
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PREDICTED TOTAL CONCENTRATION 
FIGURE 7 

Measured vs. Predicted Hourly Average GO Concentrations for Intensive Study 



background level indicator when a westerly wind prevailed, and 
the AQAM environ input for CO was reasonable. The model predicted 
that CO concentrations due to environ sources were nearly constant 
over the entire airbase. 

To check the validity of AQAM predictions for CO emissions due 
to aircraft, trailer 1 measured concentrations (now assumed to be 
reasonable background CO) were subtracted from the measured con- 
centrations at trailers 3 and 4. Figure 7 shows good agreement for 
the very limited data available. The higher predicted aircraft CO 
values at trailer 3 may result either from inaccurate specification 
of aircraft idle CO emissions in the hot refueling area or from a 
too slowly-spreading plume. It was observed at NAS, Miramar that 
many times aircraft were operated at above idle RPM in modes tradi- 
tionally input into the model as "idle". To briefly examine the 
effect of the specified engine RPM on the model predictions, all 
"idle" operations were changed to "normal". For most engines this 
decreases the relative amounts of CO and UHC while increasing NOX 
and PM. However, total emittants increase since the fuel flow rate 
increases. In reality the "normal" setting is far too high to be 
realistic. This high throttle setting is more realistic for CO 
and UHC (since they generally decrease rapidly with RPM at low speed 
and then level off) than for NOX and PT (which generally increase 
in a linear manner with RPM). A more accurate method would have bee 
to use "h idle + h normal", or etc. The increased power setting 
improved the comparison with measured CO at trailer 3 (Fig. 7). A 
change of 1 in stability category input to AOAM could also 
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significantly change the predicted concentrations at trailer 3. 

In addition to predicting reasonably accurate concentra- 
tions at specific receptors, a model should also correctly 
predict concentration profiles across the receptor grid. A 
CO concentration profile across the airbase was constructed 
(Figure 8) to illustrate the variation in predicted concentra- 
tion along the wind direction. In the two cases plotted, the 
wind was from 270° and the stability category was 3. The two 
profiles were plotted along the 8 km. y-coordinate since this 
y-coordinate most nearly passed through the trailer 1-4 loca- 
tions. Predicted and measured trailer data that were available 
were also plotted. "Trailer profiles" were sketched only to 
indicate general trends and do not necessarily represent actual 
variations. The comparison shows, as expected, that the pre- 
dicted trailer 1-4 variation had a much larger gradient than the 
8 km. profile due to closer proximity to aircraft ground opera- 
tions (taxiways, hot refueling areas, parking areas). The 
measured profiles for both 2 Aug and 6 Aug were similar to the 
predicted profiles, peaking between trailers 2 and 3. The 
higher predicted values at trailer 3 were discussed above. Also 
shown in Figure 8 is the improved "trailer profile" obtained 
when engine RPM was increased from "idle" to "normal" as 
discussed above. 

2 . NOX Emissions 

Comparison of weekend/weekday data again permitted no 
significant conclusions regarding the validity of using trailer 
1 as an indicator of background NOX. 
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X COORDINATE (KM) 

FIGURE 6 

Predicted CO Concentration Profiles Along 6 km. y-coordinate 



Figure 9 presents measured versus predicted hourly- 
average NOX concentrations for trailers 1, 2, and 4. (No 
measured data were available for trailer 3 during the ten one 
hour time periods selected for validation efforts). As pre- 
viously stated, the comparison was based upon an N0 7 conver- 
sion factor for ppm to ^gm/m -5 . Predicted concentrations from 
both aircraft sources alone and total sources are plotted to 
indicate their relative magnitudes. Predicted total concen- 
trations at trailers 1 and 4 agreed with measured concentra- 
tions within a factor of approximately three. It should be 
noted that the predicted concentrations were all very small 
and varied much less than the measured data. Also, the mea- 
sured data at trailer 2 were much greater than predicted NOX 
concentrations . 

Because of the general agreement between trailer 1 mea- 
sured and predicted concentrations, it appears that trailer 1 
again provided a good representation of background concentra- 
tions. Therefore, trailer 1 measured concentrations were 
subtracted from those measured at trailers 2 and 4 and com- 
pared to predicted aircraft NOX emissions. Again, at trailer 
2 the measured (difference) values were much greater than pre 
dieted aircraft concentrations. At trailer 4 the measured 
(difference) data agreed reasonably well with predicted air- 
carft data (both were very small). Since trailer 4 and 
trailer 1 concentrations were nearly the same for both mea- 
sured and predicted data, and only approximately one-half of 
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PREDICTED CONCENTRATION 
FIGURE 9 

Measured vs. Predicted Hourly Average NOX Concentrations For Intensive Study 



the predicted trailer 4 values were due to aircraft, trailer 
4 was probably outside most of the aircraft plumes for the 
existing wind conditions. 

Because trailer 2 was located in a near-source region 
where lateral concentration gradients were large, compari- 
sons were also made to crosswind receptor concentrations. The 
(trailer 2 - trailer l) measure( ^ concentrations were compared 
to the predicted concentrations from aircraft at special 
receptor 404 (100m crosswind/south of trailer 2). The pre- 
dicted concentrations were still much less than measured 
concentrations, indicating that the predicted concentration 
gradients around trailer 2 were not enough to significantly 
improve the comparison between predictions and measurements. 

These results indicate that the NOX emissions from air- 
craft engines specified in AQAM are too low for low power 
engine operations (idle and taxi). An alternative explana- 
tion is that the aircraft engine settings for aircraft located 
around trailer 2 (hot refueling area, taxiways , and parking 
areas) are well above idle, thus producing more NOX than 
assumed by AQAM. Increasing the engine settings from idle to 
normal as discussed above did greatly improve the comparison 
between (trailer 2 - trailer 1) predictions and measurements 

*T *7 

(from 10 to 136 /igm/m -3 vs. 190 ^gm/nr measured for 2 Aug from 
1400-1500 hrs ) . 

3 . Total Hydrocarbon (THC) Emissions 

The measured versus predicted total hour ly- averaged 
THC concentrations for trailers 1, 3 and 4 are plotted in 
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Figure 10. (No measured data were available for trailer 2). 

The conversion factor used for ppm to /<gm/m^ was based on CH^ 
and was therefore only an approximation for total hydrocarbons. 

As can be seen from the figure, predicted data were significant- 
ly lower and varied much more than measured data. Measured 
trailer 1 concentrations were approximately 1.5 times greater than 
trailer 3 concentrations. This decrease is nearly the same as 
expected for downwind dispersion from far upwind sources (i.e., 
due to changes in C in equation 3) . These results indicate that 
almost all THC was probably from environ sources. AQAM predicted 
concentrations at trailer 3 were greater than those at trailers 1 
and 4 due to aircraft ground activity. If most of the measured 
concentrations of THC are in fact due to environ sources and 
measured trailer 1 values are accurate, then either AQAM values 
for THC emittants due to environ sources are low (i.e., land-use 
factors are low) or the conversion factor to>tgm/m^ is in signifi- 
cant error. The former would also imply that the values used in 
AQAM for THC emittants from aircraft sources are too high (at 
trailer 3 downwind of the hot refueling area) . This particular 
observation could have been better clarified had measured data 
been available from trailer 2. Increasing engine RPM as discussed 
above in this case reduces THC from aircraft. When combined with 
increased environ sources of THC, this change would improve the 
agreement with measurements. 

4 . P articulate (PT) Emissions 

Figure 11 is a plot of the measured (converted bscat) 
versus predicted total hourly averaged PT concentrations'. The 
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FIGURE 10 

Measured vs. Predicted Hourly Average TfIG Concentrations for Intensive 

Study 
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FIGURE 11 

Measured vs. Predicted Hourly Average PT Concentrations for Intensive -Study 



measured data were within - 401 of the mean value. The mea- 
sured values at trailers 1 and 4 were essentially the same. 

The comparison is fairly good (within a factor of three for 
7 0 ? s of the data) at trailers 1 and 4 using the aforementioned 
conversion factor for bscat to . The model, however, 

appears to overpredict PT concentrations at trailer 2. AQAM 
predicts that most of the PT concentration is from aircraft 
sources. Therefore, if trailer 1 data are good indicators of 
background PT concentration, then AQAM has low environ source 
PT input (land-use factors, vehicle mileage, etc.) and/or high 
aircraft source PT input. Increasing engine RPM as discussed 
above increases PT and makes the comparison generally less 
favorable . 

C. CONCLUSIONS AND RECOMMENDATIONS 

Approximately 50% of the predicted levels of concentra- 
tion were found to agree with measured levels within a factor 
of two. The results also indicated that: (1) predicted CO 

concentrations agreed quite well with measured data; (2) model 
predictions xvere too loxv for NOX emissions from aircraft oper- 
ating in the idle/taxi mode; and (3) predicted THC and PT con- 
centrations were too high for aircraft operating in the idle/ 
taxi mode and/or were too low for environ sources. The latter 
appears more reasonable. Agreement between model predictions 
and measured values was significantly improved by increasing 
engine RPM settings above idle in all modes normally specified 



6 3 



at idle. 



For a reasonably complete model validation to be accom- 
plished much more measured data must be obtained during a 
specific time period of observed meteorological and operation- 
al activity. The conclusions from this intensive study were 
based on very limited data and can only be considered prelim- 
inary results. Accurate data for background levels/local air 
quality are important for determination of aircraft source con- 
tributions to total emittants. It would be most beneficial to 
obtain pollution measurements on weekends at a time when air- 
craft activity is low and meteorological conditions are very 
similar to weekday conditions. If at all possible, any additional 
intensive efforts should be conducted during a period with less 
variations in meteorology. Detailed data collection should begin 
several days before the detailed operational data are collected 
in order to ensure a more complete data set than was obtained in 
this initial effort. 

Model predictions were very sensitive to the specified 
stability category. Model validation efforts could be improved 
if stability classes could be measured and specified to half- 
integer values. 

Jet dispersion rate differences from those for elevated, 
low velocity sources and variations with orientation to the 
wind direction require further study. Plume rise of jet exhausts 
should also receive additional study. 
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APPENDIX 3 



GO AND PT CONCENTRATION PROFILES FROM AIRCRAFT SOURCES 

(SENSITIVITY STUDY) 
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10.0 



AIRCRAFT CO CONCENTRATION PROFILE (RUN NO. 1 ) 
(Scale = 20 jugm/m^ per contour) 
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AIRCRAFT r? CONCENTRATION PROFILE (RUN NO. 1 ) 
(Scale = 20 per contour) 
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AIRBASE CO CONCENTRATION PROFILE (RUN NO. 1 ) 
(Scale = 1 //gm/rr* per contour) 
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AIRBASE FT CONCENTRATION PROFILE (RUN NO. 1 ) 
(Scale = 1 />gm/:n J per contour) 
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TOTAL CO CONCENTRATION PROFILE (RUT? NO. 1 ) 
(Scale = 50 p gm/m^ per contour) 
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(Scale = 20 /Jgm/xn j per contour) 
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AIRCRAFT CO CONCENTRATION PROFILE (RUN NO. 2) 
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AIRCRAFT PT CONCENTRATION PROFILE (RUN NO. 2) 



83 




AIRCRAFT CO CONCENTRATION PROFILE (RUN NO. 3) 
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AIRCRAFT PT CONCENTRATION PROFILE (RUN NO. 3) 
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AIRCRAFT PT CONCENTRATION PROFILE (RUN NO. I4.) 
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aircraft co concentration profile (run no. 5) 
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AIRCRAFT PT CONCENTRATION PROFILE (RUN NO. 5) 
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AIRCRAFT ?T CON C3NT RAT ION PROFILE (RUN NO. 6) 
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AIRCRAFT CO CONCENTRATION PROFILE (RUN NO. 7) 
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AIRCRAFT FT CONCENTRATION PROFILE (RUN NO. 7) 
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AIRCRAFT CO CONCENTRATION PROFILE (RUN NO. 6) 
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AIRCRAFT FT CONCENTRATION PROFILE (RUN NO. 8) 
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AIRCRAFT CO CONCENTRATION PROFILE (RUN NO. 9) 
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AIRCRAFT PT CONCENTRATION PROFILE (RUN NO. 9) 
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APPENDIX C 

CO AND ?T CONCENTRATION PROFILES FROM AIRCRAFT SOURCES 

(INTENSIVE STUDY) 




AIRCRAFT CO CONCENTRATION PROFILE (1 AUG 1300-1i;00) 
INCREMENTED FROM 50.0 
(Scale = 50 ^gm/nP per contour) 
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AIRCRAFT FT CONCENTRATION PROFILE (1 AUG I3OO-U1.OO) 
INCREMENTED FROM 30. 0 
(Scale = 30 p gm/m^ per contour) 
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AIRCRAFT CO CONCENTRATION PROFILE (1 AUG lli.00-1500) 
INCREMENTED FROM 50.0 
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AIRCRAFT FT CONCENTRATION PROFILE (1 AUG UiOO-1500) 
INCREMENTED FROM 30.0 
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AIRCRAFT CO CONCENTRATION PROFILE (2 AUG 1 1^00-1500) 
INCREMENTED FROM 50 . 0 
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AIRCRAFT PT CONCENTRATION PROFILE (2 AUG 1 ij.OO-1 500) 
INCREMENTED FROM 30.0 
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AIRCRAFT CO CONCENTRATION PROFILE (2 AUG 1 500-1 600 
INCREMENTED FROM 50.0 
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AIRBASE GO CONCENTRATION PROFILE (2 AUG 1 500-1 600) 
INCREMENTED FROM 1.0 
(Scale = ( 4 . per contour) 
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AIRBASE PT CONCENTRATION PROFILE (2 AUG 1 500-1 600) 
INCREMENTED FROM 10.0 
(Scale = 5 ^gm/tn^ per contour) 
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TOTAL CO CONCENTRATION PROFILE (2 AUG 1500-1600) 
INCREMENTED FROM 100.0 
(Scale = 100 yogm/m^ per contour) 
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TOTAL PT CONCENTRATION PROFILE (2 AUG 1 500-1 600) 
INCREMENTED FROM 50.0 
(Scale = 50 per contour) 
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AIRCRAFT CO CONCENTRATION PROFILE (2 AUG 1515-1615) 



INCREMENTED FROM 50.0 
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AIRCRAFT PT CONCENTRATION PROFILE (2 AUG 1515-1615) 
INCREMENTED FROM 30. 0 
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AIRCRAFT CO CONCENTRATION PROFILE (3 AUG 1100-1200) 
INCREMENTED FROM 50.0 




AIRCRAFT PT CONCENTRATION PROFILE (3 AUG 1100-1200) 
INCREMENTED FROM 30.0 
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AIRCRAFT CO CONCENTRATION PROFILE (6 AUG lij.OO-1500) 
INCREMENTED FROM 50.0 
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AIRCRAFT PT CONCENTRATION PROFILE (6 AUG llj.OO-1 500 ) 
INCREMENTED FROM 30.0 
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AIRCRAFT CO CONCENTRATION PROFILE (6 AUG 1 500-1 600 ) 
INCREMENTED FROM 50.0 
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AIRCRAFT PT CONCENTRATION PROFILE (6 AUG 1 500-1 600) 
INCREMENTED FROM 30.0 
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AIRCRAFT CO CONCENTRATION PROFILE (6 AUG 1 51 5-1 61 5) 
INCREMENTED FROM 50.0 
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AIRCRAFT FT CONCENTRATION PROFILE (6 AUG 1515-1615) 
INCREMENTED FROM 30 . 0 
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AIRCRAFT CO CONCENTRATION PROFILE (7 AUG 1 500-1 600) 
INCREMENTED FROM 50.0 
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AIRCRAFT PT CONCENTRATION PROFILE (7 AUG 1 500-1 600) 
INCREMENTED FROM 30.0 
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